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AMENDMENTS TO THE SPECIFICATIONS 

Please amend paragraph 10005] as indicated: 

[0005] A conventional way to determine a localized spectrum is to use the method of the 
Short Time Fourier Transform (STFT). In STFT, tapered moving windows of the time 
domain signal are used to compute their Fourier spectra. The general form of this 
transform is: 

CO 

STFT r (/) = lf(t)g(t-T)e te *dt 

where f(t) is the time-domain function, g(t) is the window function, and e m \ s the 
Fourier kernel. This method has been summarized by Nawab and Qauatieri (see Nawab 
et ah, 1988. "Short-time Fourier transform" in Advanced Topics in Signal Processing, 
Lim, J., and Opperiheim, AL Eds.: Prentice Hall Signal Processing Series, 289-337.) and 
employed in practice by Partyka et aL (see Partyka et al, 1999, "Interpretational aspects 
of spectral decomposition in reservoir characterization: The Leading Edge, 18, 353-360.) 
and Marfurt et al. (see Marfurt et al, 2001 . "Narrow-band spectral analysis and thin-bed 
tuning: Geophysics, 66, 1274-1283,) The analysis window junction plays an important 
role in the STFT, where resolution issues are dependent to a large degree on window size. 
The longer the window size is in time, the better the resolution of the local spectrum in 
the frequency domain, but the worse the resolution in the time domain. The shorter the 
window size is, the better the resolution of the local spectrum in the time domain, but the 
worse the resolution in the frequency domain. 
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Please amend paragraph 10006] as indicated: 



[0007] In one method, for an observed sequence 



M 



where a M .m is the AR parameter m of the order AR process. €„ represents the noiso 
components. In terms of vectors, 

r - y T A 

where 



and T denotes the vector transpose. If each side of the above equation is pre-multiplied 





by complex conjugate vector X^ H and the expected value taken, then 



M M A 



where 






4 J 
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is the Toeplitz autocorrelation matrix, $ = E(XjX'.„} is the autocorrelation Amotion at 

lag time t. P M = [ p m , 0, • • • , 0]^ and p M =E J is the white noise power spectral 
density. 
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Please amend paragraph [0010] as indicated: 

[0010] The latest approach to precluding or diminishing the problems due to windowing 
involves the use of the wavelet analysis. Wavelet analysis is a newly established (since 
the late 1980s) field in mathematics and signal processing. Like the Fourier transform, 
the wavelet transform also convolves through a discrete summation or continuous 
integration the time function (signal) with a kernel function. A method of using wavelets 
is found in Chakraborty et aL (see Chakraborty et aL, 1995, Frequency-time 
decomposition of seismic data using wavelet-based methods: Geophysics, 60, 1906- 
1916.) as well as Xia (Xia, L., 1999, Spectral analysis of seismic data using wavelet 
transform: M.S. Thesis, University of Oklahoma.). In wavelet analysis, a wavelet is used 
as the kernel function in place of the Fourier kernel. For example, given a function/ft), 
its Fourier transform is: 

F{f)=]f{t)e tat dt (2) 

—CO 

where 0) = 2ft f and e~ ic * i$ the Fourier kernel. Its corresponding wavelet transform is: 

where 1/ is the complex conjugate of wavelet Wavelet \fr is a function that is square- 
integrable [y e 1} (5?)) having zero mean and localized in both time and frequency. 
Therefore, ^ represents a family of wavelets that satisfies the conditions: 
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where r ^0 and c ^ 0 . a is referred to as the scale of the wavelet and r is referred to 
as the translation parameter. Note that the wavelet is normalized such that = 1 . 
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Please amend paragraph [0012] as indicated: 

[0012] In MPD, a family of wavelets is defined by the foxm 

where T ^ 0 and <? ^ 0 . a is referred to as the scale, t as the translation parameter, and 
£ as the frequency modulation. Each wavelet in the family is called a time-frequency 
atom. If \t(t) is Gaussian, these atoms are called Gabor atoms. As shown by Mallat et al., 
Gabor atoms provide excellent time-frequency resolution, These atoms have 
combinations of all possible time and frequency widths and as a result constitute a 
redundant set. Once atoms are defined, a best match between the signal and these atoms 
is found by projecting the atoms onto the signal and then computing the maximum. A 
residue is then computed by subtracting from the original signal the product of the atoms 
and the cross product of the selected atom and the signal. This decomposition is 
continued until the energy of the residue falls below some threshold. This method has the 
ability to obtain a good resolution in both time and frequency for data within an 
intermediate frequency range like seismic data. 
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